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1 Introduction 



In the experiments performed at Jefferson Laboratory, MAMI, BNL, MIT on meson- 
nucleon, nucleon-hyperon, hyperon-hyperon reactions very rich data are accumulated. In 
order to describe the existing data, the coupling constants of pseudoscalar and vector mesons 
with baryons are needed. Calculation of these coupling constants within QCD, which is the 
fundamental theory of strong interactions, constitutes an important problem in studying 
the dynamics of the aforementioned reactions. At the hadronic scale QCD is nonpertur- 
bative and for this reason the calculation of the strong coupling constants of baryons with 
vector mesons becomes impossible using fundamental QCD Lagrangian. Therefore in order 
to calculate these constants some nonperturbative method is needed and for this aim we 
will use QCD sum rules method [1], which is more reliable and predictive in calculating the 
properties of hadrons. In the present work, the strong coupling constants of the decuplet- 
octet baryons with vector mesons g(DOV) are calculated in the framework of light cone 
sum rules method (LCSR) (more about this method and its applications can be found in 
[2]). Note that, the coupling constants of pseudoscalar meson octet baryon couplings, vec- 
tor meson octet baryon couplings and pseudoscalar meson decuplet baryons are all studied 
within LCSR in [3], [4] and [5], respectively. In this work we extend our previous works to 
investigate the g(DOV) in the LCSR framework. 

The plan of this article is as follows. In section 2 the strong coupling constants of 
DOV are calculated in LCSR method and relations between the above-mentioned strong 
coupling constants are obtained, when SU(3)f symmetry breaking effects are taken into 
account. Section 3 is devoted to the numerical analysis of the g (DOV). 

2 Light cone QCD sum rules for the vector meson 
with decuplet-octet baryon coupling constants 

In SU(3)f symmetry the coupling constants of all vector mesons with decuplet-octet 
baryons are described by the following interaction Lagrangian, 

C mt = ge l3k O\D mki V l m + h.c. , (1) 

where O D and V are the octet, decuplet baryons and V is the vector meson, respectively. 
Octet baryons and octet vector mesons are given by, 



O a « = 



n 



S ° "7e A J 



>° + > P+ K* 
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In this section we derive LCSR for the coupling constants of the vector meson with 
decuplet-octet baryons. In deriving these sum rules we start our analysis with the following 
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correlation function, 



K"° V = * j d '^ PX (V(Q) \T{v(x)fj,(0) } | 0) , (2) 

where V(q) is the vector meson with momentum q, D and O are the decuplet and octet 
baryons, 77^ and r\ are their interpolating currents respectively, and T is the time ordering 
operator. 

The sum rules for the correlation function (2) are derived by calculating it in terms of 
the hadrons from one side (phenomenological part), and calculating it in the limit — p 2 — > 
+00 in the deep Euclidean region in terms of quarks and gluons and then equating these 
representations using the dispersion relations. 

The phenomenological part can be obtained by inserting a complete set of intermediate 
hadronic states with the same quantum numbers as the corresponding interpolating cur- 
rents. After isolating the ground state contributions of the decuplet and octet baryons, we 
get the following result, 

U D ^ ov (p a) = (°l^l Q (^))(0(P2)V(g)|^(pi)) (D(pi) \%\ 0) (3) 
M p\ - m\ Pi - ml 

where O and D denote the octet and decuplet baryons, p\ = p 2 + q, P2, m 1 and m 2 are their 
four-momentum and mass, respectively, and • • • represent the contributions of the higher 
states and continuum. The matrix elements in Eq. (3) are defined as (see 



(0\r}\O(p 2 )) = X u(p 2 ,s) , 
(D(pi) \%\0) = \ D u^Pi,s) . (4) 

Using the Lorentz invariance, the matrix element (0(p2)V(q)\D(pi)) is parametrized in 
terms of three form factors g±, g 2 and g 3 as follows [7]: 

(0{p 2 )V(q)\D(p 1 )) = u(p 2 )^gi(qj - ej)-f 5 + g 2 [(p ■ e)q a - (P ■ g)e a ]75 

+ 93[(q ■ z)q a - g 2 ^j75}wa(pi) , (5) 

where P = (pi +p 2 )/2, and q and e are the vector meson momentum and polarization 
vectors, respectively. In the present work we assume the vector meson to be on shell, that 
is q ■ e = and q 2 = my. 

In order to obtain the correlation function from the phenomenological side, summation 
over spins of octet and decuplet baryons are performed, i.e., 



^2 M G°2, s)u(p 2 , s) = i> 2 + m 2 , 



E, \- ( \ (i . \( 7«7/3 ZPlaPip . Plalp ~ PWla\ 
u a (pi, s)u,(pi, s) = -(A + mi) ^ - — - + — j 



(6) 



Using Eqs. (3-6), in principle, one can obtain the expression for the phenomenological part 
of the correlation function. In doing so, however, the following two problems appear. The 
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first problem is related to the fact that the interpolating current 77^ couples not only to the 
J p = | + states, but also to the J P = \ states. The matrix element of the current 77, 
between vacuum and J p = \ states is determined by the following relation 

(o M 1~(pi)) = a( 7 „ - 4 ^t)«(pi. s ) • (7) 

2 

In deriving this relation the condition 7^p M = has been used. Therefore the structures 
with 7, at the end, or the structures ~ p±^ contain the contributions from the J p = \ 
states which should be removed, and we will not consider them in further analysis. 

The second problem is that not all Lorentz structures are independent. In order to 
overcome these problems we will use the ordering procedure of the Dirac matrices in such 
a way that it guarantees the independence of the Lorentz structures, as well as, it is free of 
the J P = \ contributions. In this work we will choose the 7^5^75 ordering. 

Taking into account this ordering and using Eqs. (3)-(6) for the phenomenological part 
of the correlation function we get 



[mi — {p + q) 2 ] (m 2 , — p 2 ) 
+ other structures . 



01 (mi + m 2 )#7 5 g M - #2#75(P • e)q ll + g^i^e^ 

(8) 



Obviously, the structures ^75g M , ^7s(p • £)q^ and 0l5 £ ^ do not contain contributions 
from J P = \ states, because the contributions of J p = \ states are proportional to 7, 
or (p + Note that in Eq. (8) we make the replacements p 2 = p and p\ = p + q. 

For calculation of the theoretical part of the correlation function from the QCD side the 
expressions of the interpolating currents of the decuplet and octet baryons are needed. 

The general form of the interpolating currents for the decuplet baryons can be presented 
as [6] 

77, = Ae abc [ {qfC lfl q^ g 3 c + (qf C 7 ^) q{ + (qfC^ql) g 2 c ] , (9) 

where a, b, c are the color indices, C is the charge conjugation operator. The values of A 
and the quark flavors q 1 , q 2 and q 3 for each member of the decuplet baryon are listed in 
Table 1. 

The most general form of the interpolating currents of the octet baryons are of the 
following form [6,8]: 



V 



= -\J\t abc (u aT Cs b ^ 5 d c - (^s aT Cd b ^ 5 u c + p(u aT C^ 5 s b y c - p^s aT C^ 5 d b ^u c 



rf = -rf + (s -+d) , 
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Table 1: The values of A and the quark flavors qi, qi and q% 



if = -rf' (s ->■ u) , 
rf" = -rf 1 (d ->■ s) , 
rf = —rf (u — ?■ s) , 



rf 



+ 



y|e abc \^(u aT Cd b ^s c + (u aT Cs b y 5 d c + (s aT C lb d 
2(5(u aT C^d b ^s c + f3(u aT C l5 s b y c + [s aT C lb d b ^ 



it 



(10) 



where (3 is an auxiliary parameter and j3 — — 1 case corresponds to the Ioffe current. There 
are the following relations between A and S° currents, which are shown in [9], 



2rf° (d^s)+rf° = -V3ri A , 
2rf° (u-> s)+ r] 110 = V3r] A . 



In principle, having the explicit form of the interpolating currents, the correlation func- 
tion can be calculated straightforwardly from the QCD side. Before calculating the cor- 
relation function from the QCD side we will try to find the relations among the invariant 
functions for the structures ? M , iflbip • £)Qh and ^756^. The relations between the 
invariant functions are structure independent, while their explicit terms are structure de- 
pendent. In establishing the relations between invariant functions, we will follow the works 
of [3-5]. The main power of this approach, which we present below, is that it it takes 
into account £77(3) / symmetry violating effects. We will further show that all correlation 
functions which are needed for the determination of the coupling constants of the vector 
mesons with decuplet-octet baryons can be written in terms of only one invariant function 
for each structure. 

Following the works [3]-[5], let us consider the correlation function describing the S*° — > 
S°p° transition and show that this transition is described by only one invariant function. 
The obtained result is enough to establish relations among E*° — > S°p° and S* + — > S + p° 
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and £*~ — > £~p° transitions. As has already been noted, the relations among invariant 
functions are structure independent, i.e., the relations are the same for all structures. The 
correlation function for the S*° — > £°p° transition can formally be written as 



n s* ^£ ( V = 5V>o«« n i( M > d ' s ) + 9/fldd^i(u, d, s) + g p o ss U 2 (u, d, s) . (12) 
The meson current is represented as 



u.d.s 



and for the p° meson g p o uu = —g p o dd = l/y/2 and g p o ss = 0. The functions III, and n 2 
describe the emission of the p° meson from u, d and s quarks, respectively, and are defined 
formally as follows: 

IIi(u,tZ,s) = <mm|S°S*°|0) , 

U 2 (u,d,s) = <ss|S°S*°|0) , (13) 

It follows from Eqs. (9) and (10) that, the interpolating currents of S*° and S° are 
symmetric with respect to the exchange u -H- d, and for this reason H[(u, d, s) = Hi(d, u, s). 
Using this result Eq. (12) can be written as 



V2 



U!(u,d,s) - ni (d,u,s)] . (14) 



Obviously, in the exact isospin symmetry limit n E *°~*' s0p0 = 0. 

The invariant functions responsible for S* + — > S + p° transition can be obtained by 
replacing d ->■ u in U^u, d, s) and using T,*°(d ->■ u) — ^/2Y,* + and S°(rf ->■ u) — -\^2T, + , 
which gives 

-2<mm|S + S*+|0) =4n 1 (w,w,s) . (15) 

The factor four appearing on the right-hand side of Eq. (15) is the result of the fact that 
there are four possible ways that the p° meson can be emitted from the u quarks. Using 
Eq. (15) we obtain the following relation for the invariant function III, which is responsible 
for the transition S* + — > S + p°, 

= -v^IIi(«, u,s) . (16) 

The relation for the £*~ — > £~p° transition can be obtained by making the replacement 
u ->■ d in Eq. (12) and using the fact that £*°(w ->• d) = v^X*" and £°(w ->■ d) = ~V2T,~, 
as a result of which we get, 

n E*-->E- P ° = _^2 ni (d, d, s) . (17) 
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Under exact isospin symmetry limit, from Eqs. (16) and (17) we get, 



n S* + ^S+p« = n E*-^E-p« (lg) 

We proceed now to calculate the invariant functions involving A resonances. For this 
aim we consider A + — > pp° transition. Using the fact that A + = E* + (s — > d) and p = 
— E + (s — > d), we obtain from Eq. (16) that, 

n A+ ^° = -(v««(H£ + s* + |o))( s ^d)- (g P o ss (ss\^ + i:* + \o))(s^d) 

(u,u,d) -=U 2 (u,u,d) . (19) 

v2 

Obviously, one can easily see that 

n A ° V = V2U ± {d, d, u) - -^Il 2 (d, d, u) . (20) 

v2 

Following similar lines of reasoning we obtain, 

U^^" = ^=U 2 (s,s,u), 

Tf*-^-p" = -Lll 2 (s,s,d). (21) 

The remaining relations between the correlation functions involving p, oj and mesons 
are presented in Appendix A. 

Up to this point the relations involving neutral p meson are obtained. We can now 
try to obtain similar relations among the invariant functions for the transitions involving 
charged p meson. For this purpose let us consider the matrix element (dd |E°E*°| 0), where 
d quarks from E° and E*° baryons form the final dd state, and the remaining u and s quarks 
are being the spectators. In the same manner, in the expression (ud |E + E*°| 0), d quark 
from E*° and u quark from E + form ud state and the remaining u and s quarks remain as 
spectators. Therefore, these matrix elements should be related and the explicit calculations 
indeed confirm this expectation. We find that, 

n E*<^ s+p - = ^ | S+S * | ^ = _^ (dd |E°E*°| 0) 

= -v^ni(d,u,s) . (22) 
Making the exchange u ■<-)■ d in the above expression, we obtain, 

n E*°_ s - p+ = <J M | S -E*o| ) = v^<^|E°E* |0) 

= y/m 1 {u,d,s) . (23) 

Using these arguments and performing similar calculations, we obtain the following 
relations among invariant functions involving charged p mesons: 
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(24) 

The relations among the invariant functions involving K* and mesons can be obtained 
easily using the similar arguments. These relations are presented in Appendix A for a given 
Lorentz structure. 

All of the obtained results for the coupling constants of DOV can qualitatively be 
understood from a simple diquark+ quark picture of the baryons in the following way. Let us 
consider the Hermitian conjugate channel, i.e., O + V — > D. Denote O as 0(qiq 2) (fe) where 
qiq 2 forms the diquark and q^ is a single quark. In a reaction, for example pp + — > A ++ , the 
single ^3 (in this case d) is the V-absorbing quark, and therefore, this reaction is described 
by U q3 (in our notation H2). In the reaction np + — > A + , the V-absorbing quark is from 
the diquark (here dd), and this reaction is described by, let us say, IT. 

For the reactions of A decays with the participation of p°, both functions IT and Tl 2 
contribute. Obviously, only Tlx and Yi 2 contribute to the reactions S* — > Xp° and S* — > Hp , 
respectively. Similar situation take place for the reactions involving u and <p mesons. 

From calculations the following relation between IT and H2 is obtained 

n 2 (w, d, s) = — II 1 (s, u, d) — Hi(s, d, u) . 

This relation leads to the result that all coupling constants of the vector mesons with the 
decuplet-octet baryons can be written in terms of only one invariant function without using 
the SU (3) f flavor symmetry, which is the main result of this work. 

We now concentrate on calculating the invariant function Hi. For this aim, the corre- 
lation function which describes the transition S*° — > S°p° is enough. In deep Euclidean 
region, — p\ — > 00, — pi — )■ 00, as we have already mentioned, the correlation function can 
be evaluated from the QCD side using OPE. In order to obtain the expressions of the 
correlation functions from QCD side, the propagator of the light quarks and the matrix 
elements of the nonlocal operators q(xi)Tq'(x 2 ) and q(x\)G ^q' \x 2 ) between the vacuum 
and the vector meson states are needed, where V represents the Dirac matrices relevant to 
the case under consideration, and is the gluon field strength tensor. 

Up to twist-4 accuracy, matrix elements (V(q) \q(x)Tq(0) \ 0) and (V(q) \q(x)G ^(0)1 0) 
are determined in terms of the distribution amplitudes (DA's) of the vector mesons [10-12]. 
These DA's are presented in Appendix B. 

In further analysis, we use the following expression for the light quark propagator 

S q (x) = ^-^- { ^(l-^A-—rnl(qq) (l-^) 
qK J 2tt 2 x 4 4tt 2 x 2 12 V 4 *J 192 m/ \ 6 V 



n 



n s*°^-p4 



T j£*+^Ap+ 



n 



A+- 



►np^ 



-2ni(d,s,s) , 
^2 



2Hi(u, s, d) + Hi(u, d, s) 



2Ui(u, u, d) , 
v^2rii(d, u, s) , 

-n 2 (s, s,u) , 

-2Ui(d,d,u) . 



2Ui(d, s, u) + u, s) 
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%9s I du {m^ G ^ ux)(y>lv ~ ux * g ^ ux ^a^ x 2 



32tt 2 



■x 2 A 2 
In ( — ) + 2 7E 



(25) 



where 7^ is the Euler constant, A is a scale parameter, and we will choose it as a factor- 
ization scale, i.e., A = 1.0 GeV (for more detail, see [13]). In the calculations, SU(3)f 
symmetry violation effects are included in the nonzero strange quark mass and strange 
quark condensate. These effects are also taken into account in calculation of the DA's 
[10-12]. 

Having the expressions of the light quark propagator and those of the DA's, the theo- 
retical part of the correlation functions can be calculated. Equating both representations of 
correlation function and separating coefficients of Lorentz structures ^75^, iin^{P' £ )% and 
H>lh £ ixi and applying Borel transformation to both side of the correlation functions on the 
variables p 2 and (p + q) 2 in order to suppress the contributions of the higher states and con- 
tinuum (see [14]), we get the sum rules for the corresponding vector meson decup let-octet 
baryon couplings. The contributions of higher states and the continuum are subtracted 
using quark-hadron duality. After standard calculations, for each Lorentz structure the ex- 
pressions for the invariant functions LT^ are obtained and their expressions are presented 
in Appendix C. Here superscript a refers to the invariant functions relevant to the 
coupling constants g±, #2 and #3, correspondingly. 

For a given transition D — > OV, once the Borel transformed and continuum subtracted 
coefficient functions n^, n^ 2 ^ and are obtained, the coupling constants can be written 
as 

1 1 m i I m i I m i- 

9l = e 1 ^ ^ M i+ M 2 m j , 

m 1 +m 2 \o^D 

2 2 2 

g 2 = e 1 ^ +J ^ + M 'i+ M 2 llf } , 

2 2 2 

g 3 = e ^ + ^ + M^ + Mg n (3) _ ( 26 ) 

m V ^O^D 

It follows from these expressions that for the vector meson decuplet-octet baryon strong 
coupling constants, the residues of baryons are needed. The residues of baryons are obtained 
from the analysis of two-point correlation function as are given in [6,8, 14]. The currents 
of the other baryons can be obtained from S° current by making appropriate substitutions 
of quarks. For this reason, for determination of the residues, we give the sum rule only for 
S° and 

A| oe -"4o/M* = _|1_ (5 + w + 5m(x) _ + fi)2 {Uu) {dd) 



111? 



1QM 2 



) (ss) ( (uu) + (dd) 
m s ^ (mm) + (dd) j + (m u + m d ) (ss) 
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32vr 4 " 



M 2 
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2 L 5m 2 , 

3 V 72M 2 
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rrin 



12M 2 



(27) 



where x = sq/M 2 , and 



i=0 



The contribution of the higher states and continuum in the invariant functions are 
subtracted by taking into account the following replacements 



e-</^ 2 M 2 (ln^- 7B J 
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3 Numerical analysis and discussion 



This section is devoted to the numerical analysis of the sum rules for the vector meson 
octet-decuplet baryon coupling constants. The main input parameters of the light cone 
sum rules in our case are the vector meson DA's. The DA's of the vector mesons are given 
in [10-12]. The values of the leptonic constants fy and fy, and of the twist-2 and twist-3 
parameters a|, a^,(Jy, ^ly, ^lv> 4v w 3v 4v K zvi °°ivi ^zvi as wen as twist-4 parameters 
Ci) ^4) Ci-, Ci~i K \v> K iv are gi ven i n Table (1) and Table (2), respectively, in [12]. The value 
of the other input parameters which are needed in the sum rule are (qq) = —(0.243 GeV) 3 , 
m 2 = 0.8 [8], (g 2 s G 2 ) = 0.47 GeV 4 [1]. 

It should be noted here that, the masses of initial and final baryons are close to each 
other. Therefore we can choose Mf = Mf = 2M 2 , and consequently uq = 1/2. Hence, in 
further numerical analysis, the values of the DA's only at uq = 1/2 are needed. 

It follows from the explicit expressions of the sum rules for the vector meson decuplet- 
octet baryon coupling constants that, in addition to the DA's, they also contain three 
auxiliary parameters, namely, Borel mass parameter, continuum threshold Sq, and the pa- 
rameter (3 in the interpolating current. Therefore, we need to find the working regions 
of those parameters where the results of the vector meson decuplet-octet baryon coupling 
constants are practically independent of them. 

The upper bound of the Borel parameter M 2 can be found by requiring that the higher 
states and continuum contributions to a correlation function should be less than, say 40% 
of the total value of the same correlation function. The lower bound of M 2 can be found by 
requiring that the contribution of the highest term with the power of 1/M 2 be 20-25% less 
than that of the highest power of M 2 . Using these restrictions, we obtain the working region 
for the Borel parameters. The continuum threshold is varied in the regions s — (m B + 0.5) 2 
and s = (m B + 0.7) 2 . 

In Fig. (l)-(3) we present the dependence of the couplings g±, g2 and g% for the S* + — > 
E + p° transition at five different values of the parameter /3 at sq = 4 GeV 2 . We observe 
from these figures that these couplings have good stability in the "working" region of M 2 . 
Obviously, the coupling constants are also expected to be independent of the auxiliary 
parameter /3. In order to find the working region of /3 where these couplings, we depict 
in Figs. (4)-(5) the dependence of g±, g-i and g% on cos# for the S* + — > S + p° transition, 
where is defined as tan# = /3. It follows from these figures that the common working 
region of j3 for the coupling constants is —0.5 < cos0 < 0.3, where these constants exhibit 
weak dependence on (3. Up to now we apply the standard procedure in analysis of the sum 
rules, i.e., Sq is chosen to be independent of M 2 and q 2 . However, in general the continuum 
threshold should depend on M 2 and momentum transfer squared q 2 [15]. This dependence 
leads to, so called, systematic uncertainties, and for this reason, the usual criteria for the 
stability of the results on M 2 does not provide ordering of the realistic errors. We consider 
these systematic errors to be about 15% (see [15]). Under these conditions, the result of 
our analysis predict that g 1 = 6 ± 2, g 2 = 2 ± 0.5 and g 3 = 20 ± 4. 

The same coupling constants for the other transitions are presented in Table 1 (here 
we present results only for the absolute values of the coupling constants). Only those 
results which could not be obtained by the SU (2) symmetry rotations are presented in this 
table. The errors presented in Table-1 take into account the uncertainties coming from 
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the variation of M 2 , so and (3, as well as uncertainties inherit in the input parameters and 
systematic errors. 

From Tables 1-3 we get the following results: 

• As far as the coupling constant g\ is concerned, for the channels presented in Table- 
1, we observe that there is good agreement between the predictions of the general 
current and the Ioffe current for the octet baryons. All of the results presented within 
the limit of errors satisfy the relations among the coupling constants and invariant 
functions. 

• In the case of the coupling constant g 2 , there are considerable discrepancies between 
the predictions of the general current and Ioffe current for the octet baryons for 
the £*+ E°K* + , Q- ->■ E°K*-, A+ ->■ TPK* + , A++ ->• S*° ->■ AK*°, 
A~ — > T.^K*° and H*° — > Y< + K*~ channels. The characteristic property of all these 
channels is that they all involve the K* meson. Note that, for these channels, the 
relations among the coupling constants and invariant functions are also strongly vi- 
olated. In the general current case all these discrepancies can be attributed to not 
having stable region for cos 9. In other words the results for the coupling constants of 
these transitions are not reliable. For the remaining transitions the results for both 
cases are close to each other within the error limits. 

• For the coupling constant g$ our results can be summarized as follows: The relations 
between the coupling constant g% and the invariant function IIi are strongly violated 
in the channels S*~ — > Ap~, S*° — > S°p°, as well as in the channels involving meson. 
These discrepancies, as is the case for the g 2 coupling constant, can be attributed to 
the absence of the stability region g 3 with respect to the cos#. For this reason, for 
the above-mentioned transitions, the predictions for the coupling constant g$ are not 
reliable. 

The gap between the results of the two currents can be attributed to the fact that, for 
many transitions the value j3 — — 1 lies outside the stability region of j3, causing considerable 
discrepancies between the predictions of the two currents. 

In conclusion, the strong coupling constant of the vector mesons with decuplet-octet 
baryons are studied within LQSR. It is shown that all coupling constants can be calculated 
in terms of a single universal function. 
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channel 
9i 


General current 


Ioffe current 


Result 


SU(3) f 


Result 


SU(S) f 


9i 


-5.6±1.6 


-5.7 


-7.8±0.5 


-7.8 


A -+pp~ 
9i 


9.1±2.9 


8 


12.2±0.8 


11 


9l 


-5.5±1.5 


-5.7 


-7.2±0.3 


-7.8 


9i 


-10±3 


-9.8 


-13.5±0.7 


-13.5 


_A+->-E A"*+ 

yi 


-13.1±4.1 


-11.3 


-18.4+1.4 


-15.6 




8.9+2.5 


8 


12+0.5 


11 


E*+-u>.fr* 
01 


-9.5±2.7 


-8 


-12.4+0.9 


-11 


yi 


-14±4 


-13.9 


-18.8+1.4 


-19.1 


yi 


-8.6+2.5 


-8 


-11.6+0.9 


-11 


9i 


10.8±3.0 


9.8 


14.1+0.9 


13.5 


9i 


-5.1±1.4 


-5.7 


-6.9+0.4 


-7.8 


9T 


-4.9±1.3 


-5.7 


-6.5+0.5 


-7.8 


E*+^E+<A 

9i 


8.0+2.8 


8 


11.6+0.7 


11 


9i 


7.1±1.9 


8 


10.9+0.8 


11 



Table 2: The values of the coupling constant g\ for various channels. 
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_ channel 

92 


General current 
Result SU(S) f 


Ioffe current 
Result SU(3) f 


92 

A°->pp- 
92 

92 

92 

_A+->-E A"*+ 


-2.3±1 
3.3±1.5 

-2.2±1 
-4±1.7 
-16.4±6 


-2.3 
3.2 
-2.3 
-3.9 
-4.5 


-3.6±0.6 
5.3±1.0 
-3.4±0.5 
-6.3±1.2 
-27.3±7.3 


-3.8 
5.4 
-3.8 
-6.6 
-7.6 


„E*+->.S .K"*+ 
#2 


11.1±3 
-10.8±4.8 


3.2 
-3.2 


16.9±3.5 
-17.9±4.7 


5.4 
-5.4 


#2 
y2 


-17.3±6 
-6.3±4.0 
12.4±4.8 


-5.5 
-3.2 
3.9 


-27.7±6.0 
-16.1±3.4 
19.2±3.9 


-9.4 
-5.4 
6.6 


^E*+^£+w 

52 

"*0 v 770, , 

92 
92 
92 


-2.0±0.7 
-1.9±0.6 
4.4±1.8 
4.2±2 


-2.3 
-2.3 
3.2 
3.2 


-3.2±0.6 
-3.0±0.4 
7.2±1.5 
6.7±1.1 


-3.8 
-3.8 
5.4 
5.4 



Table 3: The values of the coupling constant g 2 for various channels. 
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channel 


General current 


Ioffe current 


Result 


SU(3) f 


Result 


SU(S) f 


5 3 


31.4±3.5 


12.9 


27.4±1.8 


18.7 


A°-+pp- 
53 


-33.5±3.1 


-18.2 


-35.2±2.6 


-26.4 


53 


16.2±3.8 


12.9 


27.9±1.6 


18.7 


5 3 


36.5±8.5 


22.3 


45.7±2.6 


32.3 


_A+->-E A"*+ 
53 


22.6±2.4 


25.7 


19.6±1.9 


37.3 


„E*+->.S .K"*+ 
53 


-17.5±2.3 


-18.2 


-17.3±1.8 


-26.4 


5 3 


14.6±1.7 


18.2 


12.6±2.6 


26.4 


53 


33.5±3 


31.5 


31.4±2.4 


45.7 


53 


17.3±0.7 


18.2 


16.5±1.4 


26.4 


53 


-36.5±8.2 


-22.3 


-45.5±1.7 


-32.3 


^E*+^E+w 
53 


23.9±6.5 


12.9 


27.9±2.1 


18.7 




19.1±4.3 


12.9 


23.9±1.3 


18.7 


E* + -S>E+<* 

53 


-24.2±6.2 


-18.2 


-31.7±2.1 


-26.4 


53 


-16.2±3.8 


-18.2 


-33.2±1.7 


-26.4 



Table 4: The values of the coupling constant y 3 for various channels. 
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Appendix A 



In this appendix we present the relations among the correlation functions involving p, 
ui, 4> and K* mesons. 

• The vertices involving p meson. 



n s*°^Ap° = — 1 -[2U 1 (u,s,d) + U 1 (u,d,s)-U 1 (d,u,s)-2U 1 (d,s,u)} , 
V6 

u A-^np- = 2^U 1 {d,d,d) , 

n A++^ PP + = _2 v ^n 1 (w,w,w) . (A.l) 
The vertices involving K* meson. 



u a+^°k*+ = -^/2[u 1 (s,u,d)-U 1 (s,d,u)\ , 

u a^ak*+ = ^ [ n 1 ( a ,d,«)-n 1 ( a ,«,d)] , 

n A0 - E -** + = -2n!( s , d,d), 
nE . + _ s o*. + = 2ni(u>s >u ) ; 

n A++ - Bf ^ + = -2>/3n 1 ( s ,«,«), 

n s *°^*~ = v^n^M) , 

n n-^nox*- = _y2n!( M ,rf,s) , 

n E *"^*" = 2 V / 2n 1 (s,d,d) , 
n s *°^ s+ **- = -2^(5, d,d) , 

n s--A*r- = _^[ 2ni ( u , s ,d) + ni(M,s)] , 
n s *°- s °** = v^n 1 (d,«, s ) , 
n H *°" Ar ° = ^n^d^ + ri^d,^*)] , 

v3 

n a--.E-A-o = 2 n 1 (d,d,s) , 

n E -°-^ = - v ^n 1 (d,d,«), 

n n-^H-K*o = 2v / 3n 1 (s, s,s) , 



n E*+^*° = _ 2ni ( s , u , u ), 
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n A-->E-tf*° = _ 2v ^n i (rf,rf,rf) , 

n s*-^s-i?*o = _2U 1 (d,s,d) , 

nA o_ s o^o = _^ [Ui ( s d >u ) + ni ( S , U ,d)] , 

v3 

n A+^S+^o =2ni(s>u>u) (A2 ) 

• Finally, we present the expressions of the vertices involving u and mesons. 



n E *°^ = -^ r [U 1 (u,d,s) + U 1 (d,u,s)\ , 
v2 

n s*+^s+o, = _y2n 1 ( M)M ,s) , 
n s*-^s-u, = ^^(^^S) , 
n A+^ = _v^[ni(d,u,u)-ni(u,u,d)] , 

n A °^ = >/2[ni(u, d, d) - n^d, d, u )] , 
n s *°^ s ° w = -v^ni(«,s,s) , 

n s*°^Au, = — Lpn^w, s , d) - 2n!(d, s, u) - n^d, «, s ) + n^w, d, s )] , 

V6 

n S*"^EV = _[ ni ( s , d, u) + n^s, U, d)] , 

n E * + - E+ * = 2n 1 (s, u , u ), 
n E *^ E ^ = -2^(5, d,d) , 

■qA+— -qA — ^ra<j!) q 

n s*<W~V = 2Il 1 (s,s,u) , 

n H*-^H-<^ = _ 2 n 1 (s, s ,d) , 

n s*°-+A4> = -l = [u 1 (s,u,d)-U 1 (s,d,u)] . (A.3) 
v3 

Obviously, in the 517(2) symmetry limit, n A+ ^ AX * + , n A °^ Ai? *°, n A+ ^, U A °^ nuJ , 
■qE*°^Aw an( j -qE*°-5>A0 are a jj e q Ua i ^ Q zero Note that, in calculating H D ^ ov using Eqs. 

(A.1)-(A.3), one should use the wave functions of the corresponding vector meson in U 1 . 
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Appendix B 



In this appendix we present the DA's of the vector mesons appearing in the matrix 
elements (V(q) \q(x)Tq(0) | 0) and (V(q) \q(x)G^q(0)\ 0), up to twist-4 accuracy [10-12]: 



(V(q, A) |gi(x) 7M g 2 (0)| 0) = f v m v { q, / due™** 

q-x Jo 



e x -x 



2 2 



■ x -x \ r 1 
q-x J J 



due iUqx g v ± (u) 



1 e x -x 



— -x 



2" "(q-.r) 2 

1 ^"/^A, 



ml \ due taqx 



g 3 (u) + <l>\\{u)-2gl(u) 



(V(q,\)\qi(xh, l5 q 2 (0)\0) = --e^e x q a x p f v m v J due^g\{u) , 
(V(q,X)\q 1 (x)a,M^M = -ifiWfa - e%) f due™** 



2 2 

, Tfl\rX a / \ 



+ 



e x -x 
(q-x) 2 



(q^x v - q u x^) [ due iUqx 
Jo 



I i 

h\ - -<f> ± - -h 3 (u) 



m 



2 /•! 



+ / duetuqx hM-Mu) 



(V(q,X) \qi(x)a al3 gG^(ux)q 2 (0)\0) = / 



T 2 
V m V 



e -x 



2q-x 



q-x jo 



q a q^9fsu - mtfav - q^g^ + gpq v g ait 



+ fWv 



-A„_L 



-A„_L 



,A L 



x Jva^^a, 



+ fWv 



-A„-L 



-A_L 



A„± 



_A„_L 



q^a9pv - q^p9av - q^ a 9^ + q^^g^ 



X 



+ 



fv m2 



V 



q-x 



q a qp£pX v - q^q^e a x u - q^e^x^ + qpq v e a Xp 



+ 



x 



q-x 



q a q^ v xp - qpq^e u x a - q^e^x^ + q^e^a 
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(V(q,\) \qi(x)g s G^(ux)q 2 (0)\0) = -ifim v (efa - e%) J Va t e^ a " +ua ^S(a t ) , 
V(q,X) |gi(x)&G^(ux)75g 2 (0)|o} = -ifim v {e*q v - e%) J Va^^^S^) , 
V(q,X)\q 1 (x)g s G flu (ux) lal5 q 2 (0)\o) = f v m v q a (efa - e%) J Va ie ^ +ua ^ A(a t ) , 
(V(q, A) \Ux)g s G^(ux)i la q 2 {Q)\ 0) = f v m v q a (efa - e%) [ Va t e^ +ua ^V(a t ) , (B.l) 



where G^ u = {l/2)e iXl/a pG a ^ is the dual gluon field strength tensor, and J Vai = J da q daqda g 5(l- 
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Appendix C 



nS 1} = - 12/iWl + P) + 2rM]^(«o) 

+ flm v {\ - P)[8i 3 (A, l-v) + 12^ v (« ) - ^vW - 2B(u )]} 
+ (lE ~ In ^) { (K(l + 3/3) - 2m*/3]i 2 (S, 1) 

+ [m s (3 + 0) - 2m d ]i 2 (S, 1) + 2[m d (l - P) + m s (l + 2/3)]i 2 (7i, 1) 
+ 2[m d (l - /3) - m s (2 + P)]i 2 (T 2 , 1) + [2m d /3 + m s (l + P)]i 2 (T 3 , 1) 

+ [2m d p-m s (l + 3P)]i 2 (%,l)) 

j v m v f[2m d /3 + m s (l + /3)]ii(75, 1) + [2m d /3 - m s (l + 3^)]ii(75, 1)) 



24x/3tt 2 

- 38^^ H [2md " ms(3 + ~ 2M 2 [2m,/3 + m s (l + /3)]#(«o)) } 



my 



■{3/^m y (^ 2 G' 2 )[2m (i/ 9 - m s (l + 3/3)]ii(7a, 1) 



+ 4tt 2 



1728\/37r 2 M 4 I 

4((ss)Km s rr4(l - (5) - 6(g 2 s G 2 )m d P] + (dd)[m 2 m d m 2 v (l - P) 

- 6(g 2 G 2 )m s f3] + mgm^(l - f3)[m d {dd) + m a (ss)]) ^y(^o)} 

- 2pm d i 2 (Tu l + 2v)- 2m s (l + p)i 2 (Ti, v) + m s (l + 3p)i 2 (Ti, 1) 
+ [m s (3 + p)- 2m d ]i 2 (T 2 , 1) + 4[m d (l + p) - m s }i 2 (T 2 , v) 

+ [m s (l +p) + 2Pm d ]i 2 (T 3 , l-2v) + 4:Pm d i 2 (%, v) 

- 2m s [i 2 (%,l - v) + pi 2 (%,2 - v)}} 

- 432 ^ M2 {24tim 4 vi2P(dd) - (1 + 3P)(ss)U%, 1) 

+ Qflm 3 v (l - P)(^[m d (dd} - m s (ss}}i 2 (V,5 - Av) - [m d (dd) + m s (ss)]i 2 (A, 5 - 6vf) 

- flm 2 m v {l - p)[m d (dd) + m s (ss)][2M{u ) + ^(n )] 

+ 2j\m v (2m 2 (l + p)[m d (ss) + m s (dd}} - 9m 2 v (l - f3)(m d (dd) + m a (ss»)^v(«o)} 

-55lS^{ 16 (<^ 2 )t 2 ---( 3+ «i 

- 27r 2 m 2 [(Jrf)(5 - /3) - 2(ss)(5 + 2/3)])i 4 (Cr) + 9( 5 2 G' 2 )[2 / 5m ci + m s {\ + /3)]A T (u )} 
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M 2 ( f 

~ 96^ 1^^ V ~ 3[2/?md + ms(1 + ^M^o) + ~ ms{l + W^ 5 ' 1) 

+ 4[2m d /3 + m s (l + /3)]i 2 (<S, 1) - 16[m d (l + /3) - m a ]i 2 (5, u) 

+ 4[2m d + m s (3 + 7/3)]i 2 (7I, 1) - 16m d /3i 2 (7J, l + v)- 8m s (l + p)i 2 (Ti,v) 

- 4[2m d /3 + m s {\ + P)]i 2 (T 2 , 1) - 16[m d (l + /3) - m s }t 2 (T 2 , v) 

+ 8[2m d f3 + m s {\ + p)]i 2 (T 3 , 1 - v) + 4[2m d /3 - m s (3 + 7p)]i 2 (%, 1) 
+ 8[2m d (5 + m s (l + p)]i 2 (T 4 , v) - 3[2m d - m s (3 + / 9)]i 4 (C r )) 

- 167r 2 /y [2/3 (dd) + + /3)]<^M - flm 3 v (l ~ 5 - Qv) + ^K)]} 



= {2(1 - p)[m d (dd) + m s (ss)][3B( Mo ) - 2i 3 (A, 1 - v) + ^ v (uq)] 



+ 4(1 - (3)[m d (dd) - m s {ss)]i 3 {V, 1 - v) 

- 12([^m d + (l-(3)m s ]{ss) + [(1 - (3)m d + 4/3m s ](dd))^vM} 



fv m v 



{[2/3(dd) + (1 + P){ss)][3A T (u ) + 8i 2 (Ti,v) - 4* 2 (T 3 , l-2v)- 8i 2 (%,v)} 



72V3 1 

- 8/3[2(dd) - (ss)]i 2 (S, 1) - 8[(1 + /3)(dd) - (ss)]i 2 (S, 1 - 2v) 

+ 4[2/3(dd> - (1 + 3/3)(ss)]i 2 (71, 1) + 4[2(dd) - (3 + /3)(ss)]i 2 (T 2 , 1) 

- 16[(1 + /3)(dd) - (ss)]i 2 (75,u) + 8(1 + 2p)(ss)i 2 (%, 1) + 6[2(dd> - (3 + p)(ss)]l 4 (C T ) 

- ^f([2Hm d + (1 + P)m s }H(%, 1) + 2[2/3m d - (1 + 3/9)m»]*i(U *)) } 

+ gg/^ {3(^ 2 G 2 )[2/3m, + (1 + /3)m s ] - 40^ [2/3 <d~d> + (1 + /3)<s S >]}<^(u ) , 

(C.l) 



nS 2) = (ie - ln ^) { - ~ mmd " ms) ' l(T ' 1} ~ 2m ^ l(7 ~ 3 ' 1} 

- 2(m d - m s M%, 1)] - ^£^[ 2w * - m s (3 + /3)]7 4 (B T )} 
+ [m d (ss)(l + 3/3) - m s (dd)(l + /3)]i (* , 1)} 

-|^g{K-^)(i^)^(r,i)-^(r 4 ,i)] 

- [2m d - m s (3 + /3)] (3 + (5) [h (T, u) - 2% x {%,v)\ } 
+ |^K(^>(3 + /3) - m s (dd)(l + /3)]i (*, 1) 
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+ 14 ^^ 4 (1 - (3){^[m d (dd) - m s (ss)]M<S>, 1) + i^, 1)] 

- [m d (dd) + m a (ss)][4ii($, 1) + 4ii(tf, 1 - 2u) + 3i 4 (A)] 

+ 8[2m d {dd) - m s (ss)]ii($, u) + 8[m d (dd) + 2m s (ss)]ii($, u)| 

+ 2 t^M4 K( ^ ) + m * (s ~ s)](1 " ^ )[274(B) " + ^M] 



l39 ^_ 2 u4 {(^G 2 )[6m, - 3m s (3 + /?)] - S7r 2 m 2 [(dd)(5 -(5)- 2(ss){h + 2/3)]}* 4 (B, 

+ ^^(! - /3){4K(di>i 2 (A «) - 1)] 

- 4[m d (Jd) - 2m s (ss)]i 2 (V, v) - 4m s (ss)i 2 (V, 1) 

- 3[m d (dd) + m s (ss}} [2i 4 (B) - 2i 4 (^) + ^W]} 

- ^^{a + - 2(^)][ii(r, i) - i!(r 4 , i)] 

- [2(dd) - (ss){?> + P)\[h(T,v) -hfav) + 3l 4 (B T )]} 

- ^ 6 ^ 2 (1 - /3)[4< 2 (A l-v) + 4* 2 (V, l-v) + 2u(M) - 2i 4 (4) + ^vM] 



{2[m d (l - 3/3) + m s (l + 30)]h(T, 1) + 4[2m d /3 - m s (l + 3/3)] ?1 (T 4 , 1) 



24\/37r 2 

+ [2m d - m s (3 + P)][2h(T,v) - 4^(7^) + 3? 4 (B T )] - 4m d (l - $)h{%, 1)} 
f"m 3 r 

+ - £){ 8 [ 4 *o(*, 1 - 2u) + 4i (tt, 1) + u) - ii($, 1 - 3u) 

-i 1 (*,l-2 U )]-3i 4 (A)} , (C.2) 



nS 3) = - 0)[2z 3 (A 1 - u) - 2^ 3 (V, 1 - u) - 2? 4 (B) + 2i 4 (4) 

+ 5^y(%)] " 6/^[2^m d + (1 + /3)m s ]0^( % )} 

+ (7- - ^ ^) { " |^ (K(l + 3/3) - 2m^]< 2 (5, 1) 

- [m s (3 + P) - 2m d ][i 2 (S, 1) - i 2 (%, 1)] + 2[m d (l - /3) - m s (2 + /3)]i 2 (71, 1) 
+ 2[m d (l + m s (l + 2P)]i 2 (T2, 1) - [2m, + m s (l + /3)]i 2 (T 3 , 1)) 



+ ^0}f 2 (2<[2m d - m s (3 + /3)]i 4 (B T ) + M 2 [2m d /3 + m s (l + /3)]<^M 
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(3[m d {ss) + m s (dd}](m 2 + 2M 2 )^. v (u ) 



144 V3M e 

+ 2304^2^ {t 2m ^ + m ^ + « A rM " MTs, 1 - 2v)] 

- 8[m d (l +/3)- m s ][i 2 (S, l-2v)- i 2 (T 4 , 1 - 2v)] - 8(2m d - m s )i 2 (S, 1) 
+ 4[m s (3 + 0)- 2m d ]i 2 (Tx, 1) + 16[m d (l + /5) - m a ]i 2 (7I, u) 

- 4[2m d /3 - m s (l + 3/3)]i 2 (75, 1) - 8[2/3m, + m s (l + /3)]* 2 (T 2 , v)} 

~ 2g g^ 2 (l - «{Wil> - m.(f a >][8< 1 ($, 1) + 9i!(*, 1)] 

- 2[m d (dd) + m 8 (ss)][4ii($, 1) + 4ii(tf , 1 - 2u) + 3i 4 (A)] + 16[2m d (dd) - m s (ss)]ii($, u) 
+ lQ[m d (dd) +2m s {ss)]i 1 ($,v) + lQ[m d (dd)i 2 (A, 1) - m s (ss)i 2 (A,v)] 

- 16[m d (dd)i 2 (V,l) - m s (ss)i 2 (V, v)] j 

" 2 ^^ 2 {2(1 - ^)K(Jd) + m s ( S - S )][7 4 (B) -T 4 (4)] 

(ss)[2m d (l + f3)- m s (l - /?)] - (dd)[m d (l - 0) - 2m s (l + /9)])^ v (« )} 



fv m 



432V3tt 2 M 2 



(^G 2 ) [6m d - 3m s (3 + /?)] - 87r 2 m 2 [(dd}(5 - 0) - 2(ss)(5 + 2/3)]}i 4 (l 

ge jpm^ + m s (l + /3)][3A r («o) - 4* 2 (S, 1) - 4* 2 (T, 1) 

+ 8* 2 (T 3 , u) + 4* 2 (T 4 , 1)] + 16[m d (l + /3) - m s ] [i 2 (5, u) + i 2 (7I, u) - * 2 (T 4 , u)] 

- 4[6m,/3 - m s (l + 5/3)}i 2 (S, 1) - 8(m d + m s )(l + /3)* 2 (T 3 , 1)} 

- ^ 2 ^ 2 (1 - «) - <i($, 1 - 3t;) - <!(*, 1 - 2«) 

- i 2 (V, 1 - u) + i 2 (A 1 - v)] - 3T 4 (A)| 

+ ^^{ 3 <[ 2 ™<i - m s {3 + £)]i 4 (B T ) - 4n 2 [2(dd)P + + /3)]<^K)} 

- - ^){K(Jd)i 3 (^, 1) - m s (ss)i 3 (A «)] - K(dd)i 3 (V, 1) - m s (7ss)t 3 (V,v) 

- 3[m d (dd) + m s {ss)][u(M) - h(<f>v)]} 
f v my ([m d (l -(5) + 8m s f3]{dd) + [8(5m d + m s (l - ^)](^>}^W 



24\/37r 2 

- ^-^|[2/3(^) + (1 + /3)(ss)][4 ?2 (T 3 , 1 - 2u) - 3A T ( Mo )] + 24[2(Jrf) - (ss)(3 + /3)]i 4 (Br) 
+ 8[(dd)(l + /3) - (ss)][i 2 («S, 1 - 2u) - i 2 (T 4 , 1 - 2v)\ + 8/3[2(dd> - (ss)]i 2 (S, 1)} 
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fy™y 
18^3 



[2{dd) - (3 + p)(ss)\i 2 (Ti, 1) - 4[<dd)(l + /3) - (ss>]i 2 (7J, u) 



+ 



[2/?(Jd) - (1 + 3/9)<ss)]i 2 (75, 1) + 2[2(Jd>/3 + (1 + /3)(ss)]i 2 (75,u)} 



- {3(9 2 sG 2 }[2m d p + m,(l + - 407r 2 mg[2(Jrf)/3 + (ss}(l + /3)]}<#M ■ (C.3) 



The functions i„, i 4 and ? 4 are defined as 

io(<f>,f(v) 
ii((f>,f(v) 
«2(0,/(u) 

U(f(u) 
u(f(u) 



Von / dv(f)(aq,a q ,ag)f(v)(k-u )9(k-u ), 
Jo 

Dai \ dv4>(aq,a q ,a g )f(v)9(k-u ) , 
Jo 

J Dai J dv<j)(a q ,a q ,a g )f(v)5(k-u ) , 
J Vcti J dv(f)(a q ,a q ,a g )f(v)S'(k-u ) , 
duf(u) , 

du(u - u )f(u) , 



no 



where 



k = a q + a„f , mq 



, M\Ml 



M = 



Ml + Mf ' 



Note that d, s), (i = w, d, s) do not depend explicitly on the u quark properties 

m u and (uu). The dependence is implicit in the meson distribution amplitudes and leptonic 
constants. 
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